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Abstract 
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1— I , 

[ We study effective mass generation of off-diagonal gluons and infrared abelian 

> : 

, dominance in the maximally abelian (MA) gauge. Using the SU(2) lattice 

m ■ 

QCD, we investigate the propagator and the effective mass of the gluon field 

| in the MA gauge with the U(l)3 Landau gauge fixing. The Monte Carlo 

G\ ; 

' simulation is performed on the 12 3 x 24 lattice with 2.2 < j3 < 2.4, and 

also on the 16 4 and 20 4 lattices with 2.3 < < 2.4. In the MA gauge, the 
r~| , diagonal gluon component shows long-range propagation, and infrared 



X 



abelian dominance is found for the gluon propagator. In the MA gauge, the 



off-diagonal gluon component behaves as a massive vector boson with the 
effective mass M a g ~ 1.2 GeV in the region of r > 0.2 fm, and its propagation 
is limited within short range. We conjecture that infrared abelian dominance 
can be interpreted as infrared inactivity of the off-diagonal gluon due to its 
large mass generation induced by the MA gauge fixing. 

PACS number(s):12.38.Gc, 12.38.Aw, 11.15.Ha 
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I. INTRODUCTION 



The quark-confinement mechanism is one of the most important subjects in the nonper- 
turbative QCD (NP-QCD) and in the hadron physics. From a phenomenological point of 
view, quark confinement is characterized by the linear inter-quark potential V(r) ~ or with 
the hadronic string tension a ~ lGeV/fm, which is obtained from the Regge trajectories [JTJ] 
of hadrons and also from the lattice QCD simulation M. This confining force arises from the 
one-dimensional squeezing of the gluonic flux between the quark and the anti-quark, which 
is actually observed in the lattice QCD simulation [BJ. 

On the quark- confinement mechanism, Nambu [|J first proposed the dual-superconductor 
picture in 1974 using the concept of the electromagnetic duality [||] in the Maxwell equation 
and the analogy with the one- dimensional squeezing of the magnetic flux as the Abrikosov 
vortex in the type-II superconductor ||. In the dual-superconductor picture, color-magnetic 
monopoles are assumed to condense, and then the color-electric flux between the quark and 
the anti-quark is squeezed as a one- dimensional tube due to the dual Higgs mechanism 

HI- 

The possibility of the appearance of monopoles in QCD was pointed out by 't Hooft in 
1981 using the concept of the abelian gauge fixing ||, which is defined by the diagonalization 
of a gauge-dependent variable as G su(N c ). In fact, the abelian gauge fixing is a 

partial gauge fixing which remains the abelian gauge degrees of freedom on the maximally 
torus subgroup {/(l)^ -1 C SU(N C ). As a remarkable feature in the abelian gauge, color- 
magnetic monopoles appear as the topological objects relating to the nontrivial homotopy 
group, U 2 (SU(N C )/U(1) N ^) = Z^- 1 [9-12]. 

Then, the quark- confinement mechanism can be physically interpreted with the dual 
Meissner effect in the abelian dual Higgs theory [10,11,13-15] under the two assumptions 
of abelian dominance P JT0|Jl^JT6[| and monopole condensation [17-19] in the abelian gauge. 
Here, abelian dominance firstly named by Ezawa and Iwazaki in 1982 |Tf| means that only 
the diagonal gluon component plays the dominant role for the NP-QCD phenomena like 
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confinement. In 1990 's, the approximate relations of abelian dominance for quark confine- 
ment [18,20-22] and chiral symmetry breaking [p3| , [2~4] are numerically shown in the lattice 
QCD Monte-Carlo simulation in the maximally abelian (MA) gauge [18,20-26], which is a 
special abelian gauge in case of N c = 2. 

Abelian dominance is usually discussed on the role of the diagonal (abelian) component 
of the gluon field. However, in terms of the off-diagonal gluon, abelian dominance can be 
expressed that off-diagonal gluon components are inactive at the infrared scale of QCD 
and can be neglected for the argument of NP-QCD. Here, abelian dominance as infrared 
inactivity of off-diagonal gluons may be interpreted with a large off-diagonal gluon mass. In 
this paper, we study the gluon propagator and the effective mass of off-diagonal gluons in 
the MA gauge, using the SU(2) lattice QCD Monte-Carlo simulations. 

II. INFRARED ABELIAN DOMINANCE AND MASS-GENERATION 
HYPOTHESIS ON OFF-DIAGONAL GLUONS IN THE MA GAUGE 

In the maximally abelian (MA) gauge, the nonabelian gauge symmetry SU(iV c ) of QCD 
is reduced into the abelian gauge symmetry U(l) 0-1 , and accordingly diagonal gluons can 
be distinct from off-diagonal gluons. For instance, diagonal gluon components behave as 
neutral gauge fields like photons, and off-diagonal gluon components behave as charged 
matter fields on the residual abelian gauge symmetry PJl0|. 

In the MA gauge, abelian dominance for quark confinement and chiral-symmetry break- 
ing has been shown as the approximate relation in the lattice QCD Monte-Carlo simulations 
[18,20-24], and only diagonal gluons seem to be significant to the infrared QCD physics, 
which we call "infrared abelian dominance". In this section, we consider infrared abelian 
dominance and the off-diagonal gluon property in the MA gauge in the SU(2) QCD. In 
terms of the off-diagonal gluon, infrared abelian dominance means that off-diagonal gluon 
components are inactive and their contribution can be neglected at the infrared scale of 
QCD in the MA gauge. 
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As a possible physical interpretation for infrared abelian dominance, we conjecture that 
the effective mass of the off-diagonal gluon = -^(^ ± iAjJ is induced in the MA 
gauge. If acquires a large effective mass M Q s, the off-diagonal gluon propagation is 
limited within the short range as r ^ M~ s , since the massive particle is propagated within 
the inverse of its mass. Hence, in the mass-generation hypothesis of A^ in the MA gauge, 
the massive off-diagonal gluons cannot mediate the long-range interaction and their direct 
contribution is expected to be negligible at the infrared scale, although there remains the 
indirect contribution of off-diagonal gluons through the coupling to the diagonal gluon. 
(This is similar to the negligible contribution of heavy mesons for the nuclear force at a 
large distance, where only light mesons as pions play the significant role.) 

The mass-generation hypothesis on off-diagonal gluons in the MA gauge is formally 
expressed in QCD as follows. The QCD partition functional in the MA gauge is expressed 



MS 



Sqcd = / VA,ex V {tS QCB [A^}5($± A [AJ)A FP [A^ (1) 



where App[/1 M ] denotes the Faddeev- Popov determinant |p8 |. Here, $^ A [A^] denotes the 
off-diagonal component of $ma[^] = [D^, r 3 ]], which is diagonalized in the MA gauge 
T2| . ^8|p^| , and therefore the MA gauge fixing is provided by 5($ma[^a»])- The mass- 



generation hypothesis of off-diagonal gluons A^ is expressed as 

^qcd = / VAlexp{iS Abel [Al}} J VA±^{iS™[A±]}F[A,]. (2) 

Here, S^fA^] denotes the U( 1)3- invariant action of the off-diagonal gluon with the effective 
mass M s as 

S&[Afi = J d A x {-- (Df™A+ u - D?*A+) {D^AL - D^ hcl A^) + M 2 oS A^] . 

(3) 

with the U(l)3 covariant derivative D^ heX = <9 M + ieA 3 ^. Here, SAbei[^] is the effective 
action of the diagonal gluon component, and ^""[A^] is a U(l)3-invariant smooth functional 
in comparison with exp {^^[-A^]} at least in the infrared region. 



Since the massive off-diagonal gluons cannot carry the long-range interaction and would 
be inactive at the infrared scale, it seems natural to expect that only diagonal gluons are 
propagated over the large distance and carry the long-range confining force in QCD. In 
fact, infrared abelian dominance may be explained with the mass generation hypothesis of 
off-diagonal gluons. (This interpretation of infrared abelian dominance may be similar to 
the dominant contribution of photons for the large-scale interaction in the Weinberg-Salam 
model, where only the massless photon field is propagated over the large distance and the 
massive weak bosons does not contribute to the large-range interaction.) 

Since the effective mass M Q s is closely related to the off-diagonal gluon propagation, 
especially on the interaction range, we study the gluon propagator in the MA gauge using 
the lattice QCD Monte-Carlo simulation. In the following sections, we aim to estimate the 
effective mass of the off-diagonal gluon from the numerical result of the gluon propagator in 



the MA gauge [18,27 



III. MASSIVE VECTOR BOSON PROPAGATOR 

In this section, we try to understand the behavior of the massive vector boson propagator 
in the Euclidean metric from a simple model, for the preparation of the analysis on the 
effective gluon mass in the MA gauge in the lattice QCD. We start from the Lagrangian of 



the free massive vector field with the mass M ^ in the Proca formalism |29] 



C = l -(d,A u - d u A,) 2 + hvl^A^ (4) 
in the Euclidean metric. The propagator G^ u (k; M) of the massive vector boson is given as 

in the momentum representation H50f| . The Euclidean propagator G^ v {x — y; M) in the 
coordinate space is obtained by performing the Fourier transformation as 

G^{x — y; M) = (A,(x)A u (y)} = J ^ e ik <^G^(k; M). (6) 
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For the investigation of the interaction range, it is convenient to examine the scalar-type 
propagator G J w (r;M), since G^{r;M) depends only on the Euclidean four- dimensional 



distance r = y (x^ — y^) 2 - As shown in Appendix [A], the scalar-type propagator G w (r; M) 
can be expressed with the modified Bessel function K\(z) as 

G^r; M) = < W ^)> - / (gi^FTM^ ( 4 + w) 

^ {X - V) T^ + ^(x -y) = -^-K l{ Mr) + ±5\x - y). (7) 



(2tt) 4 k 2 + M 2 M 2 y aj An 2 r v M 2 

In the infrared region with large Mr, Eq.(|7]) reduces to 

3VM e" Mr 



C^(r;M)^ ; 5 , (8) 
using the asymptotic expansion 

K^-^e-f r(l + n) — 0) 

for large Rez. Here, the Yukawa-type damping factor e~ Mr expresses the short-range inter- 
action in the coordinate space. Then, the mass M of the vector field A^(x) is estimated 

3 

from the slope in the logarithmic plot of ~^G^^(r; M) as the function of r, 



In { -^=G w (r; M) ) - -.1 Ir ! <-onst . ( 1 j 



This approximation seems applicable for Mr > 1 from the numerical calculation [ plf . 

For the massless free vector boson, the Euclidean scalar-type propagator in the Landau 
gauge is written as 

G M = ~ (11) 



with the four-dimensional distance r = y (x M — y^) 2 - This exhibits the Coulomb-type inter- 
action in the four- dimensional coordinate space. 



IV. MAXIMALLY ABELIAN GAUGE FIXING 



In this section, we review the maximally abelian (MA) gauge, which seems the key 
concept for the study of the dual-superconductor picture from the lattice QCD [18,20-26]. 
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In the Euclidean QCD, the MA gauge is defined by minimizing the global amount of the 
off-diagonal gluon, 

R oS ^jd 4 xtr{[D„H][D„H]^ = -jd i x £ \A«\ 2 , (12) 

Q=l 



by the SU(iV c ) gauge transformation [|T^,|l^,|^,^] . Here, we have used the Cartan de- 



composition, Ap(x) = Ap(x) ■ H + YjaL\ Nc A^ L {x)E a and the covariant derivative operator 
Dp = dp+ieAp. In the MA gauge, the SU(iV c ) gauge symmetry is reduced into the U(l) Arc_1 
gauge symmetry with the global Weyl symmetry [18,22,32-34]. and the diagonal gluon com- 
ponent Ap behaves as the abelian gauge field and the off-diagonal gluon component A® 
behaves as the charged matter field in terms of the residual U(l) 7Vc_1 abelian gauge symme- 
try. In the MA gauge, off-diagonal gluon components are forced to be as small as possible, 
and then the gluon field Ap{x) mostly approaches to the abelian gauge field Ap(x) ■ H |22| . 



Since the covariant derivative operator Dp = dp + ieAp obeys the adjoint transformation 
Dp — > QDpQ^, it is convenient to describe R Q s with Dp like Eq . (|i~2"D , for the argument of 
the gauge transformation property of -R fr- Using the infinitesimal gauge transformation, 
the local gauge-fixing condition in the MA gauge is obtained as 

[H,[Dp,[Dp,H]]]=0. (13) 

In particular for the SU(2) QCD, this leads to (idp ± eAy\ A^ = 0, or equivalently, the local 
variable $ma[^W] = [Dp, [Dp, r 3 ]] is diagonalized in the MA gauge. 

In the SU(2) lattice QCD, the system is described by the link variable Up(s) = 
expliaeA^s) 2 ^-} G SU(2) with the lattice spacing a and the QCD gauge coupling constant 
e. The MA gauge is defined by maximizing 

%A = EEtr{[/,( S )r 3 [/t(,)r 3 } (14) 

s 11=1 L J 

using the SU(2) gauge transformation. In the continuum limit a — > 0, maximizing -Rma 
leads to minimizing i? Q fj in Eq. (|T2"D . Also in the lattice formalism, the MA gauge fixing 
is a partial gauge fixing on the coset space SU(2) local /(U(l)3° cal x Weyl global ), and there 



remain the U(l)3-gauge symmetry and the global Weyl symmetry. According to the partial 
gauge fixing on SU(2)/U(1) 3 , it is convenient to use the Cartan decomposition for the 
SU(2) link-variable as U„(s) = M M (s)u M (s) with u„(s) = e ie " {s)T3 G U(l) 3 and M„(s) = 
^{e^s^+el^T 2 } G SU(2)/U(1) 3 . Here, the abelian link variable u^s), which is suggested 
to play the relevant role for confinement [18,20-22] in the MA gauge, obeys the residual 
U(l)3-gauge transformation, 

Up(s) -> ^(s)w Ai (s)^ t (s + /}), (15) 

with the gauge function u(s) G U(l) 3 . 

For the study of the gluon propagator (A°(x)A^,(y)) (a,b — 1,2, 3) in the MA gauge, the 
gluon field A a Ax) itself is to be defined as the dynamical variable by removing the residual 
U(l)3-gauge degrees of freedom. In order to extract the most continuous gluon configuration 
under the constraint of the MA gauge fixing, we take the U(l) 3 Landau gauge for the U(l) 3 - 
gauge degrees of freedom remaining in the MA gauge. The U(l) 3 Landau gauge is defined 
by maximizing 

^u ( i)L = EI>K00] (is) 

s fj,=l 

using the U(l)3-gauge transformation flip] ) [18,27,34-37]. The maximizing condition for 
-Ru(i)l leads to the Landau gauge condition d^A^ = in the continuum limit. By taking 
the U(l) 3 Landau gauge, all of the abelian link variable w M (s) mostly approach to unity, and 
the most continuous gluon configuration is realized under the MA gauge fixing condition. 

V. GLUON PROPAGATORS AND OFF-DIAGONAL GLUON MASS IN THE MA 

GAUGE IN THE LATTICE QCD 

A. Gluon Propagators on Lattices 

In this sub-section, we consider the extraction of the gluon propagator from the link 
variable U^(s) in the MA gauge with the U(l)3-Landau gauge fixing. To begin with, we 
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extract the lattice gluon field A°(s) (a = 1, 2, 3) from the link variable £7 M (s) obtained in the 
lattice QCD Monte Carlo simulation. In the SU(2) lattice QCD, U^(s) relates to A M (s) as 

/ rs+afi T a\ e T a n 

U ll (s) = VexpUeJ s dx^Afa)— j = exp [iaeAfa)— j (a = 1,2, 3) (17) 

with a small lattice spacing a, the gauge coupling constant e and the pass-ordering operator 
V. Here, ^(s) corresponds to the coarse-grained field of A a (x) on the link from s to s + a/t, 
and coincides with A"(s) approximately, when a is enough small ||. In this paper, we regard 
A^(s) as the approximate gluon field on the lattice. Here, the gluon field A®(s) is described 

as 



A l = — arctan Fj5 ( 18 ) 



G S/ELi(^) 2 



by parameterizing L^s) as U^ = U^ + ir a U® (U®, Z7* 6 R; a = 1, 2, 3) , which satisfies the 
unitarity condition, (?7°) 2 + (t/^) 2 = 1. The gluon field A^s) G su(2) is thus obtained from 
the link variable U^(s) G SU(2) generated in the lattice QCD Monte Carlo simulation. For 
the simple notation, we express A^(s) by A^s) hereafter. 

Next, we study the Euclidean gluon propagator G ab {x — y) = (A a ^(x)A b v {ij)) in the 
MA gauge with the U(l) 3 -Landau gauge. In this gauge, U^{s) is determined without the 
ambiguity on the local gauge transformation. In this paper, we study the Euclidean scalar- 
type propagator of the diagonal (neutral) gluon as 

G^\r) = (Al(x)Al(y)}, (19) 

and that of the off-diagonal (charged) gluon as 

G°Z(r) = {A+(x)A-(y)) = \ {(Al(x)Al(y)) + (A 2 (zMj (</)>} (20) 

with the off-diagonal gluons A^(x) = 4= {^4* (a;) ± 2v4 2 (x) j. These scalar-type propagators 
are expressed as the function of the four-dimensional Euclidean distance r = yj (x M — y M ) 2 . 
In G^(r), the imaginary part — | ^(Aj l (x)Af 1 (y)) — (A'j l (x)Aj l (y))^ disappears automatically 
due to the symmetry as {A^{x) A^iy)) = (A 2 (x)A 1 ^ (y)). 



B. Lattice QCD Results for Gluon Propagators 



Using the SU(2) lattice QCD, we calculate the gluon propagator in the MA gauge with 
the U(l)3 Landau gauge fixing. The Monte Carlo simulation is performed on the 12 3 x 24 
lattice with 2.2 < (3 < 2.4, and also on the 16 4 and 20 4 lattices with 2.3 < (3 < 2.4. All 
measurements are done every 500 sweeps after a thermalization of 10,000 sweeps using the 
heat-bath algorism. We prepare 100 gauge configurations at each j3 for the calculation. 

The MA gauge fixing and the U(l)3 Landau gauge fixing are performed by the iteration of 
the local maximization of -Rma and i?u(i)L> respectively, using the overrelaxation algorithm 
with the overrelaxation parameter u = 1.7 |3q| . As the accuracy of the maximizing condition 
of -Rma and -Ru(i)L, they are required to converge to 10 -7 per a sweep. 

The physical scale unit is determined so as to reproduce the string tension as o = 0.89 
GeV/fm ~ (420 MeV) 2 based on Refs. [|39H43|. In Fig.0 and Table 1 ||-||, we summarize 
the SU(2) lattice QCD data for y/aa, the square root of the string tension in the lattice 
unit, as the function of (3 in the case of the pure gauge standard action. In the region of 
2.2 < (3 < 2.4, the lattice QCD data in Figjl] seem to lie on a straight line, and therefore 
we determine the values of \faa and the scale unit a at intermediate values of (3 G (2.2, 2.4) 
using the interior division. 

The Euclidean scalar-type propagators G^ cl (r) and G°^(r) are measured as the two 
point functions of the gluon fields on (x, 0) and (x,r) using Eqs. fll9D and (^0|). While 
the temporal direction is fixed along the longest side in the 12 3 x 24 lattice, we measure 
the correlation in all four direction for the symmetric lattices (16 4 and 20 4 ) without fixing 
the temporal direction. We show in Fig|| the lattice QCD data for the diagonal gluon 
propagator G^ cl (r) and the off-diagonal gluon propagator G°^(r) in the MA gauge with 
the U(l)3 Landau gauge fixing. Here, error is estimated with the jackknife analysis |30 



In the MA gauge, G^ el (r) and G°^(r) manifestly differ . The diagonal-gluon propagator 
G^ cl (r) takes a large value even at the long distance. In fact, the diagonal gluon A 3 ^ in 
the MA gauge propagates over the long distance. On the other hand, the off-diagonal gluon 
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propagator G°®(r) rapidly decreases and is negligible for r > 0.4 fm in comparison with 
G^ cl (r). Then, the off-diagonal gluon seems to propagate only within the short range 
as r ^ 0.4 fm. Thus, 'infrared abelian dominance' for the gluon propagator is found in the 
MA gauge. 

Finally in this subsection, we briefly discuss the finite size effect on the gluon propagators, 
G^ cl and As shown in Fig.|], the shape of G^ el slightly depends on the lattice 

parameters (/3 and the lattice size), although /^-dependence of G^ cl at (3 G [2.3,2.4] seems 
to decrease for the 20 4 lattice. (See also Figs. 3 and 4.) Then, one should be careful about the 
finite-size effect for the detailed argument on the diagonal gluon propagator Gr^ el . On the 
other hand, the single curve approximately fits almost all the lattice data of obtained 
from different lattices on (3 and the size. Therefore, the finite-size effect seems to be small 
on the off-diagonal gluon propagator in the MA gauge. (See also Figs. 3-5.) 

C. Estimation of the Off-diagonal Gluon Mass in the MA Gauge 

In this subsection, we examine the hypothesis on effective mass generation of off-diagonal 
gluons in the MA gauge. To this end, we investigate the logarithm plot of r 3 ^ 2 G°J^(r) as the 



function of r, since the massive vector boson with the mass M behaves like G^(r j ~ — ^372 



cxp(— Mr) 



for Mr > 1 [^TJ as was discussed in Section 3. 

In Fig. 3, we show the logarithmic plot of ^ 3 ^ 2 G°^(r) and r 3 ' 2 G^ el (r) as the function of 
the distance r in the MA gauge with the U(l)3 Landau gauge fixing. As shown in Fig. 3, the 
logarithmic plot of r 3 ^ 2 G°J^(r) seems to decrease linearly in the region of 0.2 fm < r < 1 fm. 
In fact, G°^(r) behaves as the Yukawa- type function 

s exp(-M off r) 
G w» ^2 ( 21 > 

for 0.2 fm < r < 1 fm in the MA gauge. From the linear slope on r 3 / 2 Gr^(r) in Fig. 3, the 
effective off-diagonal gluon mass M g is roughly estimated as 

M oS ~ 5 fm" 1 ~ 1 GeV. (22) 
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For more accurate estimation of the off-diagonal gluon mass M a g, we analyze the lattice 
data of r^ 2 G°^{r) individually at each lattice with the same (3 and the lattice size, to 
separate the systematic error relating to the lattice parameters. We show the lattice data 
for r 3 / 2 G°^(r) in Figs. 4 (a),(b) and (c) corresponding to the lattice size, 12 3 x 24, 16 4 and 
20 4 , respectively. Using the least squares method, we perform the linear fitting analysis for 
the logarithmic plot of r 3 / 2 G*(r) with the fitting range of 0.2 fm < r < 1 fm. 

We summarize in Table 2 the effective off-diagonal gluon mass M oS obtained from the 
slope analysis at each lattice (the lattice size and (3) as well as x 2 and the degrees of freedom, 
Ndf. The error of M Q s is estimated with the jackknife analysis. Here, the two crude-lattice 
data with 12 3 x 24 and (3 = 2.2, 2.23 show extremely large values on x 2 /^/, which means 
the wrong fitting of these two data, and hence the errors listed are not creditable for the two 
crude-lattice data. Except for these two wrong-fitting data, the lattice QCD data for the 
off-diagonal gluon mass in the MA gauge distribute around M oS ~ 1.2 GeV within about 
0.1 GeV error. 

We show in Fig. 5 the off-diagonal gluon mass M oS as the function of 1/V, the inverse 
volume in the physical unit, to investigate the finite-size effect on M oS . In Fig. 5, we have 
dropped the two wrong-fitting data on the crude lattice with 12 3 x 24 and (3 = 2.2, 2.23. 
The systematic errors including the finite-size effect seem to be small for M g. 

Finally in this section, we discuss the relation between infrared abelian dominance and 
the off-diagonal gluon mass M off ~ 1.2 GeV, which is considered to be induced by the MA 
gauge fixing. Due to the large effective mass M g, the off-diagonal gluon propagation is 
restricted within about ~ 0.2 fm in the MA gauge. Therefore, at the infrared scale as 
r 3> 0.2 fm, the off-diagonal gluons cannot mediate the long-range force, and only the 
diagonal gluon A 3 can mediate the long-range interaction in the MA gauge. In fact, in the 
MA gauge, the off-diagonal gluon is expected to be inactive due to its large mass M oS in the 
infrared region in comparison with the diagonal gluon. Then, we conjecture that effective 
mass generation of the off-diagonal gluon in the MA gauge would be an essence of abelian 
dominance in the infrared region. 
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VI. SUMMARY AND CONCLUDING REMARKS 



In this paper, we have studied the gluon propagator and the off-diagonal gluon mass in 
the MA gauge with the U(l)3 Landau gauge fixing using the SU(2) lattice QCD. The Monte 
Carlo simulation is performed on the 12 3 x 24 lattice with 2.2 < (3 < 2.4, and also on the 
16 4 and 20 4 lattices with 2.3 < (3 < 2.4. In the MA gauge, we have measured the Euclidean 
scalar-type propagators G^{r) of the diagonal and the off-diagonal gluons, and have found 
infrared abelian dominance for the gluon propagator. In the MA gauge, we have found that 
the off-diagonal gluon behaves as a massive vector boson with the effective mass M q r ~ 1 
GeV for r > 0.2 fm. At each lattice, we have estimated the effective off-diagonal gluon mass 
M g from the linear fitting analysis of the logarithmic plot of r 3//2 G°^(r). The off-diagonal 
gluon mass is estimated as M Q g ~ 1.2 GeV with the error about 0.1 GeV. We have also 
checked the smallness of the finite-size effect for the off-diagonal gluon mass M g in the MA 
gauge, but still now, we have to consider the finite lattice-spacing effect and the Gribov copy 



effects [p4| 446[| relating to the gauge fixing for more quantitative arguments. 

From the behavior of the diagonal-gluon propagator G^ cl (r), the diagonal gluon seems 
to behave as a light or massless particle. However, for the detailed argument on G^ cl (r), 
one should consider the finite size effect more carefully, because the diagonal gluon may 
propagate over the long distance beyond the lattice size. On the other hand, the off-diagonal 
gluon propagation is limited within the short distance as r ^ M~ s ~ 0.2 fm, and hence the 
finite size effect for is expected to be small enough, when the lattice size is taken to be 
much larger than ~ 0.2 fm. Such a tendency of the small finite-size effect on G°J^ has 
been checked using the lattice QCD simulation with the various volume. 

Due to the large off-diagonal gluon mass as M g ~ 1.2 GeV, the off-diagonal gluon 
cannot mediate the interaction over the large distance asr > -^off ; anc ^ suc h an infrared 
inactivity of the off-diagonal gluon may lead infrared abelian dominance in the MA gauge. 
Then, an essence of infrared abelian dominance would be physically explained as generation 
of the off-diagonal gluon mass M oS induced by the MA gauge fixing. As an interesting 
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conjecture, this off-diagonal gluon mass generation may provide general abelian dominance 
for the long-distance physics in QCD in the MA gauge, similar to the infrared inactivity of 
the massive weak-bosons in the Weinberg-Salam model. 

Finally, we consider the criterion of the scale where abelian dominance would hold in 
the MA gauge. We conjecture that ~ 0.2 fm would be regarded as the critical scale 
on abelian dominance in the MA gauge. For the short distance as r < M^ 1 ~ 0.2 fm, 
the effect of off-diagonal gluons becomes significant, and hence all the gluon components 
have to be considered even in the MA gauge. On the other hand, at the long distance as 
r ^> M^ 1 ~ 0.2 fm, the off-diagonal gluon cannot mediate the interaction and would be 
negligible in comparison with the diagonal gluon, which may lead abelian dominance in the 
MA gauge. In this way, abelian dominance is expected to hold for the long-range physics as 
r ^> M~ s , and abelian dominance would break at the short distance as r < AQ 1 in QCD in 
the MA gauge. 
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APPENDIX A: FORMULA ON EUCLIDEAN PROPAGATOR 

In this appendix, relating to Eqs.(|7|) and (pT|), we add several formula on the vector boson 
propagator in the Euclidean metric. For the massive vector case, the scalar-type propagator 
G^(r;M) = (A ft (x)A IM (y)) is expressed as 

G tv (r; M) = 3/ 0-^^ + ±- 2 5\ X - y). (Al, 
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Since G^(r; M) depends only on the four- dimensional distance r = J (x M — y^) 2 , one takes 
Xfi — Vfj, = {r, 0,0,0) without loss of generality for the calculation of G w (r; M). Then, the 
integration in Eq. (|Al|) can be expressed with K\(z) as 



d 4 k I r d 3 k ( r°° dk ikor 1 



f Jk-jx-y) 1 = f 

J (2txY k 2 + M 2 J 



-e 



(2tt) 4 ~ k 2 + M 2 J (2tt) 3 \J-oo 2tt " k 2 + k 2 + M 2 J 



2 „ 

Vfc 2 +M 2 

27r ) 3 2Vfe 2 + M^ ~4tt 2 ./o """VWTm*' 

— / dEe~ E WE 2 - M 2 = — -M 2 / cfor^v 7 ^! 
47r 2 Jm Ait 2 Ji 

IJiK^Mv) (A2) 
47r^ r 



with = y k 2 + M 2 and e = E/M. Here, we have used the integration formula for the 
modified Bessel function, 

/oo 
dte~ z Wt 2 - 1 > 0). (A3) 

Thus, the scalar-type propagator G M(U (r; M) can be expressed as 

GW(r; M) = ^y^i(Mr) + ^5 4 (x - y). (A4) 

For the massless vector case, we use the covariant Lagrangian in the Euclidean metric, 

£ = \ (d,A v - d v A,f + ^(Vy 2 (A5) 

with the gauge-fixing parameter £. The massless vector-boson propagator is obtained as 

^w^ilv+ce-i)^} (A6) 



p i y fc 2 

in the momentum representation [3D|. The scalar-type propagator G^x) in the coordinate 



space is written as 

3 + £ 1 



47r 2 r 2 

which becomes Eq.([TTl) for the Landau gauge £ = 



= ( A7 ) 
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FIGURES 

FIG. 1. The square root of the string tension in the lattice unit, \/aa, as the function of (3 in 
the SU(2) lattice QCD with the pure-gauge standard action, using the data in Refs. [38-42]. In the 
region of 2.2 < (3 < 2.4, the lattice QCD data seem to lie on a straight line. 

FIG. 2. The SU(2) lattice QCD results for the scalar-type gluon propagators 

G^ c \r) = (Al(x)A*(y)) and G^(r) = (A+(x)A-(y)) as the function of r = ^{x, - y^ in 
the MA gauge with the U(l)3 Landau gauge fixing in the physical unit. The Monte Carlo simula- 
tion is performed on the 12 3 x 24 lattice with 2.2 < (3 < 2.4, and also on the 16 4 and 20 4 lattices 
with 2.3 < (3 < 2.4. The diagonal-gluon propagator G^ el (r) takes a large value even at the long 
distance. On the other hand, the off-diagonal gluon propagator G°^{r) rapidly decreases and is 
negligible for r > 0.4 fm in comparison with G^ el (r). 

FIG. 3. The logarithmic plot of r^G^Jr) and r^G^ c \r) as the function of the distance r 
in the MA gauge with the U(l)3 Landau gauge fixing, using the SU(2) lattice QCD with 12 3 x 24 
(2.2 < (3 < 2.4), 16 4 and 20 4 (2.3 < (3 < 2.4). The off-diagonal gluon propagator behaves as the 
Yukawa-type function G°£(r) ~ cxp( ~^ offr) with M oS ~ 1 GeV for r > 0.2 fm. Therefore, the 
off-diagonal gluon seems to have a large mass M r ~ 1 GeV in the MA gauge. 

FIG. 4. The logarithmic plot of r§G°^(r) and r§G^ cl (r) in the SU(2) lattice QCD with (a) 
12 3 x 24 and 2.2 < (3 < 2.4, (b) 16 4 and 2.3 < (3 < 2.4, (c) 20 4 and 2.3 < (3 < 2.4 in the MA gauge 
with the U(l)3 Landau gauge fixing. The slope corresponds to the effective mass, and the effective 
off-diagonal gluon mass is estimated as M ff ~ 1.2 GeV for r ^ 0.2 fm both on these lattices. The 
finite lattice-size effect seems to be small on the off-diagonal gluon propagator. 



20 



FIG. 5. The off-diagonal gluon mass M Q g as the function of 1/V, the inverse lattice volume 
in the physical unit. The off-diagonal gluon mass M Q g is obtained from the slope analysis on 
the logarithmic plot of ?"2G*(r) at each lattices ((3 and the lattice size) with the fitting range of 
0.2 fm < r < 1 fm. In this figure, we drop the two wrong-fitting data on 12 3 x 24 with (3 = 2.2 
and 2.23, because these data show extremely large values on \ 2 (See Table 2). 
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TABLES 



Table I. The SU(2) lattice QCD data for y/aa, the square root of the string tension in 
the lattice unit, at various j3 for the pure-gauge standard action. These data are taken from 
Refs. [38-42]. N s and N T denote the lattice size iV s 3 x N T . We add the physical length of a 
to reproduce y / a=0.89GeV/fm. 

Table 2. The effective off-diagonal gluon mass M Q s obtained from the slope analysis on 
the logarithmic plot of r^G°^(r) at each lattice {(3 and the lattice size) with the fitting range 
of 0.2 fm < r < 1.0 fm. Here, \ 2 an d the degrees of freedom, N d f, are also listed. Except 
for the first two wrong-fitting data with large x 2 , one finds 1.1 GeV ^ M Q s < 1.3 GeV. 
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